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Acoustic Source Modeling for High Speed Air Jets 


M.E. Goldstein and Abbas Khavaran 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


The far field acoustic spectra at 90° to the downstream axis of some typical high speed jets are calculated from 
two different forms of Lilley’s equation combined with some recent measurements of the relevant turbulent source 
function. These measurements, which were limited to a single point in a low Mach number flow, were extended to 
other conditions with the aid of a highly developed RANS calculation. The results are compared with experimental 
data over a range of Mach numbers. Both forms of the analogy lead to predictions that are in excellent agreement 
with the experimental data at subsonic Mach numbers. The agreement is also fairly good at supersonic speeds, but 
the data appears to be slightly contaminated by shock-associated noise in this case. 


I. Introduction 

The acoustic analogy introduced by Lighthill 1 over 50 years ago remains the principal tool for predicting the 
noise from high speed air jets. Its most general formulation amounts to rearranging the Navier-Stokes equations into 
a form that separates out the linear terms and associates them with propagation effects that can then be determined 
as part of the solution. The non-linear terms are treated as “known” source functions to be determined by modeling 
and, in more recent approaches, with some or all of the model parameters being determined from a steady RANS 
calculation. The “base” flow (about which the linearization is carried out) is usually assumed to be parallel and the 
resulting equation is usually referred to as a Lilley’s 2 equation. 

The major drawback with these approaches is that the unsteady effects, which actually generate the sound, must 
be included as part of the model. This clearly puts severe demands on the modeling aspects of the prediction, which 
usually amount to assuming a functional form for the two-point time-delayed velocity correlation spectra. These 
predictions should, however, be less sensitive to the details of the model when it is possible to neglect variations in 
retarded time across the source correlation volume. It is therefore fortunate that this seems to be a reasonable 
approximation when performed in an appropriate moving frame of reference, 3 assuming, of course, that the Mach 
number is not too large. The source models are usually tested by comparing them with measurements of the far field 
acoustic spectrum at 90° to the downstream jet axis, which is believed to be uninfluenced by propagation effects. A 
major purpose of this paper is to show that this spectrum can be accurately predicted by using an appropriate 
acoustic analogy approach combined with some measurements of the source function that were recently carried out 
by Harper-Bourne. 4 


II. The Acoustic Analogy Equation and its Far -Field Solution 

Reference 5 shows that the Navier-Stokes equations can be rewritten (for an ideal gas) as the Navier-Stokes 
equations linearized about a fictitious “base” flow but with different (in general non-linear) dependent variables, 
with the heat flux vector replaced by a generalized enthalpy flux and with the viscous stresses replaced by a 
generalized Reynolds stress. This is a true acoustic analogy (in the Lighthill 1 sense) in that it shows that there is an 
exact analogy bebA’een the flow fluctuations in any real flow and the linear fluctuations about a fictitious “base 
flow” due to an externally imposed stress tensor and energy flux vector. 

When the “base” flow is taken to be the unidirectional transversely sheared mean flow 


v t = 8 n U (x 2 ,x 3 ), p = p (x 2 ,x 3 ),p = p = constant 


( 1 ) 
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where x = {xj , X 2 , X 3 ] is a Cartesian coordinate system, v = { V, , V, , V 3 ] denotes the velocity, p the pressure and p 
the density, the general equations reduce to the modified Lilley’s 1 equation 


where 
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is the variable-density Pridmore-Brown 17 operator. 


c 2 =Yp/p(x 2 ,x 3 ) 


(4) 


is the square of the mean-flow sound speed, and 


Dt 3 1 3xj 


denotes the convective derivative based on U. The symbol t denotes the time, y denotes the specific heat ratio, 


/ / Y 1 / / 

Pe=P + — P' 7 ,' 7 , 


(6) 


is a generalized pressure fluctuation, 


y-1 


/ / , i o /2 , 

= -pv f v.- + -4— 5„pv + c 


i j 


(7) 
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is the generalized stress tensor, and 


<= ~pv'h' 0 - q' + Cj v'j 


( 8 ) 


is the generalized stagnation enthalpy flux. 
Here, 


v' = Vi-b n U, 


h>h' + l -v\ 


(9) 


( 10 ) 


and Ji denote fluctuating quantities, with h being the enthalpy 0 (/ being the viscous stress, and and q f being 

the fluctuating viscous stress and heat flux vector respectively, which are believed to play a negligible role in the 
sound generation process and are therefore neglected in the following. 

/ 

It is usually argued that the stagnation enthalpy flux TV- , (which corresponds to the isentropic part of the 

pressure density source in the Lighthill approach 1 ) is only important for hot jets 2 ' 6 " 8 except, perhaps, at small angles 
to the downstream jet axis. 9 It is therefore neglected in the present analysis. The resulting source function not only 


significant effect on the predicted acoustic spectrum. The present research was initially undertaken in order to 
determine the appropriateness of including such a term. 

This generalized Lilley’s equation was solved in reference 10 in terms of the free space Green’s function 18 
G J,TJ , which satisfies 


du , 

contains the usual quadrupole type term, but also the dipole-like term (y 

r)x, 


D 


1 \ / 

1J which could have a 


Dt 


LG(x,t\ j,T) = 8 (x-j) 8 (f-T), 


( 11 ) 


and has outgoing wave behavior at infinity, to obtain the following expression 


P 2 {x,t 0 )= J \\y ijkl (x\y^ + iU c T 0 f o +T o )R™ l (y;Z > ,T o )dy d% dr 0 (12) 

— oo JG 


for the pressure autocovariance 20 (notice that p e reduces to p in the far field) 


P 2 {x,t 0 ) 



(x,t)p' e (x,t + t 0 )dt. 


(13) 
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The symbol V denotes integration over all space; T denotes some large but finite time interval, 


u c t q 


( 14 ) 


denotes a moving frame coordinate system, 


Yijki +T o) - J Y,y ( X, y, t } + 1 0 + i 0 ) lid (* | + n, ) dt\ , (15) 


and the propagation factor y kl (x|j,^is defined in reference 10. R^ (l r 0 ) is a moving frame correlation 

tensor, which is defined in terms of the fixed frame density-weighted, fourth-order, two-point, and time -delayed 
fluctuating velocity correlation 


1 T 

R ijki (>’• % ) = — J p v'v'j ( y , r)pv' k v' {y + x\,z + z 0 )dz 

11 _ T 


(16) 


and the second order fixed frame density weighted correlation 


1 T 

R ij ( j; ) = — J Vp v 'i (y > *) VK {y+w+\)dx 


(17) 


by 


R m {y&r 0 ) = R ijkl (M + i U C T 0 , T a )-Rjj (y;0,0)R kl (j + ^ + i U c x 0 ; 0, o) . 


(18) 


The indicated arguments refer to all three terms preceding the parentheses. 
Our interest is in the far field spectrum 



— oo 


which can be calculated by taking the Fourier transform of equation (12) and using the convolution theorem 21 to 
obtain 
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7 a>(*W = 27t J | r y (x\y;^)T* kl (x\y + ^ + iU c z 0 ;(£>)e mr °Rf kl {y,^T 0 )d^dr 0 , ( 20 ) 

— OO V 


where 


T ij j eHt ~%( x \ y,t-i)d{t-x) 

— OO 


( 21 ) 


is the Fourier transform of y.. (we use capital letters to denote Fourier transform of the corresponding lower case 
quantity) and we introduced / 0) (.Y | v) , the acoustic spectrum at x due to a unit volume of turbulence aty. i.e., 


V 


( 22 ) 


in order to simplify the formulas. The relevant far field expansion of F ( y is given in reference 10. The only 

approximation made up to this point is the neglect of the enthalpy and viscous source terms, but equation (20) will 
depend on the turbulent source correlations only through 

ftyk! {y ,\ ) s \ R w (y> (23) 

V 


if variations in retarded time across the correlation volume are neglected, i.e., if 


r li(x\y+V+iU c v,a^ 


is 


assumed to be constant over the correlation volume. 3 However, the definitions (14) and (18) imply that the 
integration variable in equation (23) can be changed back to T| , which means that 



Rij (y; 0, 0) R kl (y + tj; 0, 0 )]rf0 


(24) 


Equation (20) can now be rearranged into the simpler form 


7 (o(*M 


— > 


2k ^ 

. x J 


— sin0r..(x|j 1 )r* / (x|j ± )O*. 1/ (j,(l-M c cos0)co) ,asx^°°, (25) 
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where 



— oo 


is the spectral tensor of the source correlation and 


M 


c 



(27) 


is the convective Mach number of the turbulence. This result shows that it is only necessary to model the overall 
spectral tensor itself and not the detailed two-point time -delayed correlations of the turbulence. However, the 
radiated sound should still be relatively insensitive to the detailed turbulence structure even when the latter 
quantities are modeled (as is at least partially done below). This would not be the case if the moving frame had not 
been introduced before neglecting the retarded time variations. 3 

Our interest here is in the spectrum at 90° to the jet axis where cos0 = 0. Reference 10 shows that 


4 ( x \y) = 


(4roc) 2 |_ 


x i x j 


y-U | ((y-Q 8 i, du 
2 ij CO dy, 


x k x i V-lg 


kl 


'■(Y-l)o dU 

°\k 

co dv. 


(28) 


for 0=71/2, 


when C 


o 



= constant , i.e., in the isothermal case. 


III. The Quasi-Normal and Axisymmetric Turbulence Approximations 

To proceed further, we need to know something about the source spectral tensor O i/H The usual approach’ 12 ' 13 
is to begin by assuming that the turbulence is quasi-normal 16 (see ref. 10) in order to obtain some relations among its 
components. It then follows that 


R m {y A T 0 ) - Rij ( y ; 0 , 0) R kl {y + tj; 0 , 0) 


Rik{y^^ 0 )Rji{y^^o) 




(29) 


To further reduce the number of independent components it is usual to assume some kinematically possible 
symmetric form for the second order correlations. Early studies 23 assumed the turbulence to be isotropic, but that 
turns out to be incompatible with the Harper-Bourne 17 measurements that will be introduced below. The simplest 
assumption compatible with his results is the one introduced in references 12 and 13, namely that the turbulence is 
axisymmetric which implies that 16 
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R ,j {y,Wo) = Mftj + + c ( A£u + Do (My + M/ ) 


(30) 


where the symbols A n , B 0 , C 0 , and D n denote functions of y, T 0 , and T)j_ = yjvfe + f|? ; A 0 , B 0 and C 0 denote even 

functions of the latter quantity whil eD 0 denotes an odd function. This model is chosen because it is the most general 
of those whose mathematical properties have been studied in the literature and because it reflects the fact that the 
cross flow velocity components tend to be much more similar to one another than to the stream-wise component 
(even for non-axysymmetric flows). 

Inserting equation (30) into equation (29) and inserting the result into equation (25) via equations (24) and (26) 
yields (after a straightforward but tedious calculation that follows along the lines of the one in appendix A of ref. 12) 


4(*ly)( 4 ™) 2 


= 2 


( A 4 

CO 


0i-(y-l)0 2 + 


f y-O 2 
v 2 j 




+ 


(7-1)— |W| 


4> 


4’ 


(31) 


where 


^ J e im ° I R 22 (y> ) dl \ dx o 


(32a) 


<*> 2 J e- im °\(R^+R? 2 +Rt 2 )dx\dT 0 , 


(32b) 


<*>3 J ^{(44+24 +4+24)^ dX a , 


(32c) 


and 




(32d) 


are seemingly independent spectral functions. However, the coefficients A, B, C, and D are not all independent and, 
when compressibility effects are neglected (i.e., when p is treated as a constant), these turbulence correlations can be 
expressed in terms of two independent scalar functions of y, T 0 , T) j_ i and ty , say a and Z>, 14 ' 15 24 which scale like 
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b = pu^ B (fj / 2 


(33) 


and 


S = a ~\ \ =P u 2 D [^1’\) 


(34) 


where 


fli=r|i/A 

(35) 

fl±=rii/^i 

(36) 


L x and L l denote characteristic stream-wise and transverse length scales of the turbulence, B and G are 0(1) 
functions of the indicated arguments. 

Turbulence measurements suggest that 


8 


Al 

4Zi 


(37) 


ought to be small. In fact, Harper-Bourne’s 17 measurements (to be discussed below) suggest that 8-2.7x10 2 
Reference 10 shows that 


1/3/ 1 
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0^E~ j terms are neglected. The ratio r is defined by 


r — 



( 39 ) 


Lacking any specific data to the contrary, it seems reasonable to suppose that 



— oo 


J R 22 (>’- l l- T o )d r \d T o 

V 


= <rr) 2 <t » 0 s (rp ) 2 - f e m - f R 2 {y, i|, T„ >ft| dr c 

n J 

— oo y 


(40) 


where T is a constant. Equation (31) then becomes 


( CO ^ 


,(x|y)(27Dcc 00 ) 2 =C 0 2 O o (j,®) — co 2 +(k|VC/|) 


V'oo j 


(41) 


where 



(y-l) + 2 





(42) 


is a constant, i.e., independent of co,and 



(43) 


IV. The Harper-Bourne Spectrum 

The results cannot be made more explicit without inputting more specific information about the turbulence 
structure. This is accomplished with the aid of some recent measurements 17 of the two-point fourth-order stream- 
wise velocity correlation spectra along the centerline of the mixing layer in a low Mach number jet, which would 
most closely correspond to 


H 0 {y,% to) = - j e~' on ° (y, t], x 0 ) dx 0 

K J 


(44) 
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with the quasi-normal approximation that is being used in the present analysis. 

Harper-Bourne 17 divided H 0 into the three components (see his eqs. (2.5) and (2.7) on p. 2) 


H 0 =H 0 (y,0,o>)R 


h l± j 


( 45 ) 


where /, , l L are the spectral stream- wise and transverse length scales (not necessarily the same as the time domain 
length scales L x and L_ introduced above) and 


X 


P 


u c 


(46) 


No assumption is made about the decomposition of the correlations into products of their space and time 
components with this approach. 

The first factor can be evaluated from his measurements of ^im (y 5 0,X 0 ) , which are reasonably well 

represented by the exponential e ^ . But reference 19 shows that i? ml (y,0, T 0 ) does not in reality have a sharp 

cusp at t 0 = 0 . A better representation would therefore be 


tit (y 5 o, x 0 ) 


A/Cfi 


erfc. 


^0 

2(3 


+e >a °crfc. 




2p 


J 


2erfc.(3 


(47) 


I 

which behaves like e ' °' for large T 0 and reduces to this quantity when (3 = 0, but smoothes out the cusp at 
T 0 = 0 . It therefore follows that 20 


# 0 (j,0,co) 


4 -(pco/k)- 


Xpu x e 
7i(cxp.(3 2 crfc.[3)(A 2 +co 2 ) 


(48) 


Inserting these into equation (46) and using the result in equations (41) and (42) shows that 


4( x ly) = c ( 


Xpu[ e 


V(p®/kf 


or + 


( k | vc /|) 2 2 u%llR(y,k) 


x 2 c* (exp .(3 2 erfc. (3)(X 2 +co 2 )co 


(49) 


where 
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(50) 


R{yJ ] ) = 2n \ 

— oo 


0 


-2kit\^ /j 




and /, and l L are defined by 


T 

2nU c 


(51) 


and 


T _ 0)/ i 

1 2nU c 


Harper-Bourne obtains the best fit to his data with the non-separable form 


(52) 


R = e ~^ 


(53) 


which can be inserted into equation (50) to obtain 


R = 



(54) 


Harper-Bourne’s figure 13 shows that while \ and l L are constant at relatively low frequencies, it is the scaled 
length scales /, and l L that become constant as CO — > . The data is well represented by the functions 


2nUj 2 v ' 


(55) 


and 


= 0.15(l- e ~ 0 5St ) (56) 

2nUj v ' 


where 
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and D j denotes the jet diameter. 


V. Extension of the Harper-Bourne Data 

Unfortunately, all of Harper-Bourne’s data are taken at a single point in a very low Mach number jet, while 
practical interest is in much higher Mach number flows and the acoustic predictions require information about the 
turbulence over the entire noise producing region of the jet. We therefore attempt to extend his data by using some 
modeling assumptions along with the Wind code developed by NASA Glenn Research Center and the U.S.A.F. 
Arnold Engineering Development Center, which is a RANS code with a standard k—e turbulence model. To this 

end, we first assume that the time scale that appears in equation (49) is proportional to the k — 8 time scale y , 

i.e., we put 


r 1 



(58) 


where C x is an adjustable constant. 


Equations (55) and (56) are extended by assuming that the time and velocity scales j and U j are 


U, 


proportional to the k — e time and length scales k/ and k' 2 respectively to obtain 


^~c’k 112 

In 


1-exp 


C 


s OJk ' 

2ne 


(59) 


and 


^ = 0.3 C l k xn 
2k 


1 -exp\ -0.25 C S — 
1 27te 


(60) 


where C 1 and C S are constants. Figure 1 shows the RANS (WIND code) solution for a 2 in. diameter cold jet at the 
jet exit Mach number of 0.18 for which Harper-Bourne carried out his measurements. The two constants C' and 
c s are determined by requiring that equations (55) and (59) be in reasonable agreement with the data at the 
measurement point (x^/ Dj = 4.0, Xo / Dj = 0.50). A reasonably good approximation is obtained by selecting 
l s 

C ~ 3.0 andC ~ 0.50 . We neglected the difference in exponents between (59) and (60) but plan to remedy 
this oversight in the near future. 

The velocity ratio f is typically close to 'A in. most cases and we shall use this value in our computations. The 
constant C 0 is related to the ratio Tr (defined implicitly by eq. (40)) by eq. (43), which for y = 1 .4 becomes 

C o 2 =^(rF) 2 +0.01 (61) 
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Unfortunately, Harper- Bourne only measured the stream-wise and not the transverse velocity correlations so that V 
is essentially unknown. But if we assume, as in reference 10, that B / 2 = D = a (t|j , T 0 ) e " then this quantity 
is equal to V3 and it follows that 


Cl =0.1472. 


(62) 


We treat C x ,(3, and C 0 as adjustable constants, whose determination is described in the next section. It is necessary 
to specify the square root in equation (43) in order to fix K . Again, Harper- Bourne does not provide enough data to 
ascertain this quantity, but it becomes equal to l A in. and equation (43), therefore, becomes 


K = 



(63) 


when, as before, it is assumed that y = 1 .4 and B / 2 = D = a (ijj , T 0 ) e 


-hi/2 


VI. Comparison with Measurements 

The far-field spectra at 90° to the jet axis were calculated from equations (49) and (54) at jy = lOODj for Mach 
0.50, 0.90, and 1.5 cold jets. As in the JeNo code, 15 ' 21 the local result (49) was summed over all source points within 
the noise producing region of the jet in order to predict the actual sound field. RANS solutions for the relevant 
nozzles were obtained from the WIND code with upstream boundary conditions specified in terms of the stagnation 
pressure and temperature at the nozzle plenum entrance. The predicted turbulent kinetic energy distribution and the 
corresponding time-scale for the three jets are shown in figure 2. 

The calculated 90° acoustic spectra are compared with the subsonic SHJAR data recently acquired at NASA 
Glenn Research Center and correctly expanded supersonic data obtained at Langley Research Center in figures 3 
through 5. Atmospheric attenuation was removed from all measurements in order to make a lossless comparison 

with predictions. We expect the parameters C l and C S , which were determined from the Mach 0.18 jet RANS 
solution, to be independent of both Reynolds number and Mach number. The adjustable constants C x , (3, and C Q 
were determined by obtaining the best fit with the Mach 0.5 data. The resulting values of C T and (3 turn out to be 

0.35 and 0.10 respectively. The scale factorC 0 turns out to beCj- 0.133, which is in remarkably good 
agreement with the value calculated from equation ( 61 ). 

The calculated spectra are in excellent agreement with the subsonic data over the entire frequency range. The 
agreement is not quite as good for the supersonic case, but it is likely that this data contains a small amount of shock 
associated noise that is not accounted for by the theory. 

Equation (63) implies that k = 0.27 when C Q -0.133 . Figure 4 (with K = 0.9 ) shows that there is very little 
difference between the results obtained with K given by this equation and those obtained with K = 0. Equation (49) 
shows that I ~ CO 4 as CQ -A 0 in the latter case, which is consistent with the conventional wisdom that, at least 
for cold jets, the sideline noise is dominated by a quadrupole- type source as originally proposed by Lighthill. 


VII. Discussion 

A simpler and, we believe, more elegant form of Lilley’s equation was derived in reference 19 by introducing 
the new dependent variable 
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to obtain 


where 


L ^IL*L^JL d A 

Dt dxj dxj 3.T[ 


OX i OX; 


(65) 


( 66 ) 


L is defined by equation (3) and v' is defined by equation (9). In this context it is usual to neglect the dipole-like 
, dn 

dx ; 


term (y-1 )h'^~ rather than the dipole-like term ^y_ for cold air jets. It was also shown in 


Dt 2 dx, 


reference 19 that this dipole-like term will not even appear in equation (66) if c 2 is replaced by c 2 +(c 2 ) (where 
|c 2 j = (y- l)// denotes the fluctuation in the squared sound speed) in the operator L defined by equation (3), so its 


neglect can also be interpreted to mean that the sound speed fluctuations have a negligible effect (relative to the 
mean) on the acoustic propagation. Since equations (2) and (65) are both exact, any differences in the predictions 
must be attributable to the neglect of these terms. 

When the preceding analysis is applied to the present equations (i.e., equations (64) to (66) rather than 
equation (2)), the final result is still given by equation (41) but withK = 0, a slightly different definition of the 
density weighted source correlations, and 


Ct 



(67) 


which becomes 


C 2 ~ 0.187, 


(68) 


when the values T and ^ obtained in the previous section are inserted-a result that is fairly close to the previous 
value. The principal difference between these predictions is therefore due to the second term in the factor 


to 2 + ( k | V £/|) 2 


, which does not significantly affect the high frequency behavior of the solution but has the 


potential of causing I m to exhibit the dipole-like behavior 


I m ~ to 2 as co — > 0 


(69a) 


in the formulation discussed in this paper rather than the quadrupole-like behavior 
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as co — > 0 


(69b) 


CO 


that always occurs when the predictions are based on equations (65) and (66). But the computations and data 
comparisons of the previous section show that the second term in this factor is relatively small for cold jets and that 
good agreement is achieved independently of whether that term is included. This is because the low frequency roll 
off of the acoustic spectrum is primarily determined by the peak frequency distribution of the local spectra and not 
by their low frequency asymptotes. We note, however, that the value of C 0 given by equation (62) is slightly closer 
to the “fitted” value than the one given by equation (68), but-given the uncertainty of the approximations used in the 
source modeling-this difference is not large enough to distinguish between these two forms of the acoustic analogy. 


VII. Concluding Remarks 

The research was initially motivated by the need to distinguish between the two forms of the acoustic analogy 
described above. Unfortunately the results turned out to be inconclusive-with both forms of the analogy yielding 
excellent agreement with the data. Our hope is that similar comparisons for hot jets or jets with more complex flow 
fields will provide the required selectivity. But until this is done, our recommendation would be to base the jet noise 
predictions on the formulation (65) and (66), as was done in reference 21, since this leads to much simpler formulas- 
especially at angles other than 90°. 

Finally, it is worth noting that the adjustable constant (3, which measures the curvature of the temporal 
autocovariance at x = 0 , is relatively small and is therefore consistent with experimentally observed turbulence 
spectra. -0 It is, however, somewhat puzzling that the high frequency roll off of the predicted acoustic spectra turns 
out to be fairly sensitive to this parameter. It is also rather unfortunate, because this quantity is difficult to measure 
with any accuracy. 
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Figure 1. — Predicted turbulent kinetic energy (top), and frequency scale 
(bottom) in a Mach 0.18 cold jet. r = | 
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Figure 2a. — Predicted turbulent kinetic energy (top), and frequency scale 
(bottom) for a 2 in. diameter cold jet at Mach 0.50. r = |jjJ . 


Figure 2b. — Predicted turbulent kinetic energy (top), and frequency scale (bottom) 
for a 2 in. diameter cold jet at Mach 0.90. r = I jj_ I . 
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Figure 2c. — Predicted turbulent kinetic energy (top), and frequency scale (bottom) 
for a Mach 1.50 convergent-divergent nozzle 
with 1.68 in. exit diameter, r = |j_J . 



Figure 3. — Spectrum at 90° and at R/D 3 = 100 for a Mach 0.50 cold jet. 
Prediction (dashed line); data (solid line), / = , R = |j| . 
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Figure 4. — Spectrum at 90° and at R/D,= 100 for a Mach 0.90 cold jet. 
Prediction with K — 0.0 (dashed line); K — 0.90 (dash-dot); 
data (solid line), / = ^ , R = |^| . 



Figure 5. — Spectrum at 90° and at RJDj= 100 for Mach 1.5 cold jet. 
Prediction (dash-dot); data (solid line), / = , R = | j| . 
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